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Abstract: The dynamical responses of random held Ising model at zero tempera¬ 
ture, driven by standing magnetic held wave, is studied by Monte Carlo simulation 
in two dimensions. The three diherent kinds of distribution of quenched random 
held are used here, uniform, bimodal and Gaussian. In all cases, three distinct dy¬ 
namical phases were observed, namely, the pinned, structured and random. In the 
pinned phase no spin hip is observed. In the structured phase standing spin wave 
modes are observed. The random phase is shown with no observed regular pattern. 
For a hxed value of the amplitude of the standing magnetic held wave, in the region 
of small quenched held, the system remains in a pinned phase. In the intermediate 
range of values of random held, a standing spin wave mode (structured phase) is 
observed. The regular pattern of this spin wave mode disappears for higher values 
of random held yielding a random phase. The comprehensive phase baundaries are 
drawn in all three cases. The boundary of pinned phase are analytically calculated 
for uniform and bimodal types of quenched random helds. 
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1. Introduction: 

The random field Ising model (RFIM) is a simple model to understand various 
physical phenomena, like, hysteresis[l], Barkhausen noise [2] observed in ferromag¬ 
netic materials, avalanche[3], return point memory effects etc. The RFIM was stud¬ 
ied theoretically with intense attention. RFIM was solved exactly[4] on Bethe lattice 
and the qualitative behaviour of magnetisation as a function of external magnetic 
field was found to depend on the coordination number of Bethe lattice. The effect 
of coordination number on the nonequilibrium critical point in RFIM was studied[5] 
in details. Hysteresis in RFIM with asymmetric distribution of quenched random 
fields in the limit of low disorder was studied[6] and the spin flip was found to be 
related to bootstrap percolation. Athermal hysteresis was also studied[7] in antifer¬ 
romagnetic RFIM recently. The statistics[8] and the dynamical critical behaviour 
of avalanches[9] was studied in RFIM. All these studies mentioned here are static in 
a sense that the time dependence of any quantity is studied. 

The dynamical behaviours of the driven RFIM are also quite interesting. The 
RFIM has also served to study the dynamics of domain wall. The motion of domain 
wall by magnetic field in a 2D ultrathin Pt/Co/Pt film (showing perpendicular 
anisotropy and quenched disorder which is analogous to RFIM) was studied[10] 
experimentally by magneto optical polar Kerr imagning technique. The motion 
of moving interface was analyzed numerically[ll] in the RFIM driven by external 
magnetic field. Later, the depinning transition of driven interface was studied[12]. 
The creep motion of interface in driven RFIM was studied [13] and the nonlinear 
field-velocity relationship was found. 

All the studies mentioned above have a very common feature. Firstly, the ex¬ 
ternal driving magnetic field was varied is a quasi-static manner. Secondly, the 
driving magnetic field is uniform over the lattice at any given instant of time. The 
RFIM also shows interesting dynamical behaviours if the uniform driving field varies 
rapidly. The dynamic hysteresis and dynamical athermal phase transition was al¬ 
ready studied[14] in RFIM. More interesting news comes from the fact when RFIM 
is driven by a field having both spatio-temporal variation. Very recently, the RFIM 
was studied[15] in the presence of propagating magnetic field wave. Various nonequi¬ 
librium phases are observed[15] and a phase boundary was drawn. In this paper, the 
dynamical responses of RFIM driven by standing magnetic field wave is studied in 
details. The paper is organised as follows: the model and simulation are discussed 
in section-2, numerical results are reported in section-3 and the paper ends with a 
summary of the work in section-4. 


2. Model and simulation: 

The Ising model on a square lattice {Lx L), in the presence of quenched random 
magnetic field and driven by a magnetic field having spatio-temporal variation, can 
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be described by the Hamiltonian, 

H = - JSs(a;, y, t)s{x', y\ t) - y)s{x, y, t) - T.hs{x, y, t)s{x, y, t). (1) 

Where s{x, y, t) = ±1 is Ising spin at site (x,y) and at time t. J > 0 is ferromagnetic 
spin-spin interaction (nearest neighbour only) strength, hr{x,y) is quenched random 
field and hs{x,y,t) is external magnetic field having spatio-temporal variation. In 
the present study, hs{x,y,t) represent the value of the external magnetic held at 
lattice site (x,y) and at time t. The standing magnetic wave held hs{x,y,t) can be 
represented as 


hs{x,y,t) = hoSin(27r/t)Sin(-^). (2) 

Where, ho and / represent the amplitude and frequency of the standing magnetic 
held wave respectively. Here, n represents the number of antinodes in the standing 
magnetic wave. The standing magnetic held wave spreads along the y direction only. 
The preiodic boundary conditions are applied in both directions. 

The distributions of the quenched random held hr{x,y), are considered here of 
folllowing three types: 

(a) The uniformly distributed between —w and +w, with prbability PuiK) = ^ 
if —w < hr < +w and 0 otherwise. 

(b) Having a bimodal probability distribution 

Pb{hr) = 0.55{hr — tc) -|- 0.55(h,, + w). (3) 


where 6 reprensents the Dirac delta function. 

(c) The normally distributed with 0 mean and standard deviation a {=2w here) 
with probability 


Pn{hr) 




( 4 ) 


This RFIM driven by standing magnetic held wave is studied by Monte Carlo sim¬ 
ulation at zero temperature. The zero temperature Metropolis dynamics is cosidered 
here in the manner that the spin only hips if it lowers the energy[16]. Starting from 
an initial conhguration as 


s{x,y,t = 0) = +1 V x,y. (5) 

The parallel updating rule is used here. such updates of spins constitutes the unit 
time step, i.e. Monte Carlo Step (MCS). Throughout the simulation the frequency 
of the standing magnetic held wave is kept constant (/ = 0.01). So, time period 
r = 100 MCS is required to have a complete temporal cycle of the standing magnetic 
held wave. For n loop standing wave, one would get n number of antinodes in the 
length L. As a result the wavelength of the standing magnetic held wave A = 2L/n. 
Here, the simulation is done and the results are shown for n = 4 and n = 2. The 
linear size of the square lattice is chosen L = 100 here. The wavevector, k = ^. 
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3. Numerical Results: 

Here, the dynamical responses of the model are studied for three different kinds 
of the distributions of quenched random field. For uniformly distributed (eqn-3) 
quenched random field, depending on the values ofw, ho mainly three different kinds 
of dynamical or nonequilibrium phases are observed. For ho = 1.0 and w = 1.0, the 
RFIM remains pinned, i.e., no spin flips was observed. This is called pinned phase. 
For the higher values of ho and tc, namely ho = 2.0 and w = 3.0, a structured phase 
was observed. Here, the band like sturcture of the clusters of the spins are formed. 
These bands oscillates periodically as the helds changes sign in the standing magnetic 
held wave. A typical phase of this kind may be found in a video presentation in 
http://youtu.be/5pzW3chzYnw. Note that the value of the n is taken here equals 
to 4. That corresponds to the presence of 4 antinodes in the standing magnetic held 
wave. For ho = 1.0 and w = 5.0, a random spin conhguration (with no specihc 
structure) is found. This is called random phase. Figure-1 shows snapshots of such 
dynamical phases. 

For a particular value of ho, the three phases (mentioned above) are found by 
varying the strength of the quenched random held w. The system remains in a 
pinned phase untill a value of w is reached. After that, it remains in the structured 
phase and ultimately it reaches the random phase for much higher values of w. 

All these nonequilibrium phases can be characterise by the following quatitites: 
the dynamic structure factor S{k,t) = (1/T) m{y,t)e''^^dy, where the line mag¬ 

netisation along y direction is dehned as m{y,t) = (1/T) s{x,y,t)dx. Here, the 
integrals are essentially the summation over discrete lattice point. The time average 
magnetisation over the full cycle of the standing magnetic held wave is dehned as 
Q = -^ §m{y,t)dt. Here also, the integral is essentially the discrete sum over the 
MCS. All these quantities are obtained by averaging over 100 diherent samples of 
quenched random held. 

In the pinned phase, the structure factor S{k, t) will be zero and Q will be unity. 
However, in the structured phase S{k,t) is nonzero and Q will be less that unity 
(due to the presence of spins of opposite sign). In the random phase, where no 
specihc pattern is present and approximately 50 per cent of the total number of 
spins have opposite sign. As a result, Q and S{k,t) will be nearly equal to zero. 
For a hxed value of ho the transition from pinned to random via the structured 
phase was observed by studying the S{k,t) and Q as functions of w. In the case 
of uniformly distributed quenched random held, Figure-2 shows such variations for 
ho = 2.0 and ho = 3.0. Here, the n is chosen equal to 4. For small values of w, the 
system remains in pinned phase {\S{k, y) \ = 0.0 and Q = 1.0 in hg.2(a) and hg2.(b)). 
Above a particular value of w, the system transits to the sructured ( middle of Fig-1) 
phase, where |5'(h,t)| shows a large nonzero value and Q becomes less than unity. 
This transition point can be marked by a sharp maximum and minimum of the 
derivatives and ^ respectively (hg.2(c) and hg.2(d)). This transition point 

is Wp. As the w increases, both Q and |S'(/c,t)| decreases and eventually vanishes in 
the random phase. It is observed that structured phase persists as long as |S'(/c,t)| 
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remains larger than 0.1. The structure completely disappears for |S'(A;,t)| < 0.1 
and the system gets into random phase. The transition from structured to random 
phases are marked by |S'(A;,t)| < 0.1. This transition point is called Wg- It is noted 
that, the transition values of the queched random held, Wp and Wg depend on the 
value of ho- The Wp decreases as ho increases. Physically, larger amount of standing 
wave held is required to hip the spin for weak quenched random held. On the other 
hand, Wg increases as w increases upto a certain value of w. 

The comprehensive phase boundaries (pinned-structured and stuctured-random) 
can be plotted in the plane described by ho and w. It is shown in hg-3. It may be 
noted the the bounadry of pinned and structured phases is linear. This can be 
realised as follows: the spin hip will not be possible until the local held is negative. 
The local held F = 4 —(tc + ho) remains positive as long as w + ho < 4, preventing the 
spin hip. So, the boundary of pinned phase is nothing but a straight line w + ho = 4. 
The phase boundaries are also obtained for n = 2. No considerable diherence was 
observed. 

The similar studies were done for bimodal distribution (eqn-4) of quenched ran¬ 
dom held. The qualitative behaiours are similar to that observed in the case of 
uniformly distributed quenched random held. Figure-4 shows the plots of |S'(A;,t)|, 
Q and the derivatives and ^ as functions of w for two diherent values of ho. 

The comprehensive phase baundaries (for n = 4 and n = 2) are shown in hgure-5. 
Here, the boundary of pinned-structured phase is linear and same argument may be 
applied to analyse it. However, the boundary of structured-random phase for higher 
values of ho is a line w = 4.0. As the distribution is bimodal, w can be +4 or -4. 
So, for larger values of ho, tc = 4.0 will randomise the spin structure. 

The normally distributed (eqn-5) quenched random held was also used to study 
the dynamical or nonequilibrium behavious of RFIM at T = 0. Here also, qualita¬ 
tively similar kind of behaviours were observed. Figure-6 shows the plots of \S{k,t)\, 
Q and the derivatives and ^ as functions of w for two diherent values of 

ho- However, the region of pinned phase is much smaller than that for other two 
(uniform and bimodal) cases. The quality of the data is not quite conclusive to get 
any analytic form of the phase boundary. The comprehensive phase boundaries (for 
n = 4 and n = 2) are shown in hgure-7. 


4. Summary: 

The nonequilibrium behaviours of the random held Ising ferromagnet, at zero 
temperature, driven by standing magnetic held wave are studied by Monte Carlo 
simulation in two dimensions. The three diherent kinds, namely, uniform, bimodal 
and Gaussain, of distributions of the quenched random helds are considered here. 
The random held Ising ferromagnet is driven by an additional standing magnetic held 
wave of well dehned amplitude, frequency and wavelength. The wavelength depends 
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on the number of antinodes are present in the standing wave. In the present study, 
only 2-loop and 4-loop standing waves are considered. Starting from an initial spin 
conhguration having all spins are up, the zero temperature Metropolis dynamics are 
employed. In this case, the spin only flips if its lowers the energy. The rule parallel 
updating is used here. 

In the dynamical steady state, the system shows three distinct nonequilibrium 
phases depending on the values of the amplitude of standing magnetic held wave 
and the strength of quenched random held. For a hxed value of the amplitude of 
the standing magnetic held wave, in the limit of weak quenched random held, the 
system remains in a pinned phase. In this phase, all spins remain up as considered in 
the initial conhguration. In the intermediate values of the quenched random held, a 
standing spin wave mode is observed. Here, the bands of up and down spins formed 
alternately and oscillates. This is diherent from the propagation mode studied [15] 
earlier in RFIM driven by propagating magnetic held wave. Where the bands of 
alternate up and down spins moved in the direction of the propagating magnetic 
held wave. 

The comprehensive phase boundaries are drawn in the plane formed by the 
strength of random held and the amplitude of standing magnetic held wave. These 
phase boundaries are qualitatively similar in three diherent kinds of distributions 
of quenched random held. The boundaries of pinned phases in the cases of uniform 
and bimodal distribution can be obtained analytically. 

The compound LiHoxYi_xF 4 can be modelled by random held Ising model. On 
applying standing magnetic held wave, a possible way to see these ehects experi¬ 
mentally, may be the relaxation experiments using SQUID magnetometer [17]. 


Acknowledgements: Author would also like to thank Ritaban Chatterjee for his 
help to prepare the video of standing spin wave modes. 
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Fig-1. The snapshots (taken after t = 100 MCS) of spin conhguration for uniform 
distribution of quenched random field. The black dots represent s{x,y,t) = +1. 
The pinned phase is shown (for ho = 1.0, w = 1.0) in the top. The structured (for 
ho = 2.0, w = 3.0) phase is shown in the middle. The bottom one shows the random 
(for ho = 1.0, w = 5.0) phase. Here, L = 100 and 4-loop standing magnetic field 
wave is considered in each case. The dynamics of the structured phase can be found 
in http://youtu.be/5pzW3chzYnw. 



















































































































































































Fig-2. The \S{k,t)\, Q, and ^ are plotted against w for uniform distribu¬ 

tion of quenched random held for two different values of ho and 4-loop of standing 
magnetic held wave. Diherent symbols represent diherent values of ho- (•) repre¬ 
sents ho = 2.0 and (o) represents ho = 3.0. Solid and dotted lines are just connecting 
the data points. 
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Fig-3. The phase diagram for uniform distribution of quenched random held. The 
boundaries of pinned-structured (+) phase and the structured-random (x) phase are 
shown for n = 2. The boundaries of pinned-structured (*) phases and the structured- 
random (□) phases are also shown for n = 4. The typical spin conhgurations (for 
60x60 lattice) are show in the insets for three different phases. The pinned phase 
(ho = 1.0, w = 1.0), the structured phase (ho = 2.0, w = 3.0) and the random phase 
(ho = 1.0, w = 5.0) are shown for n = 4. 
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Fig-4. The \S{k,t)\, Q, and ^ are plotted against w for bimodal distribu¬ 

tion of quenched random field for two different values of ho and 4-loop of standing 
magnetic held wave. Different symbols represent different values of ho- (•) repre¬ 
sents ho = 3.0 and (o) represents ho = 2.0. Solid and dotted lines are just connecting 
the data points. 
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Fig- 5. The phase diagram for bimodal distribution of quenched random held. The 
boundaries of pinned-structured (+) phase and the structured-random (x) phase are 
shown for n = 2. The boundaries of pinned-structured (*) phases and the structured- 
random (□) phases are also shown for n = 4. The typical spin conhgurations (for 
60x60 lattice) are show in the insets for three different phases. The pinned phase 
(ho = 1.0, w = 1.0), the structured phase (ho = 2.4, w = 2.0) and the random phase 
(ho = 1.0, w = 4.0) are shown for n = 4. 
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Fig-6. The \S{k,t)\, Q, and ^ are plotted against w for Gaussian distri¬ 

bution of quenched random held for two different values of ho and 4-loop of standing 
magnetic held wave. Diherent symbols represent diherent values of ho. (•) repre¬ 
sents ho = 2.2 and (o) represents ho = 1.0. Solid and dotted lines are just connecting 
the data points. 
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Fig-7. The phase diagram for Gaussian distribution of quenched random held. The 
boundaries of pinned-structured (+) phase and the structured-random (x) phase are 
shown for n = 2. The boundaries of pinned-structured (*) phases and the structured- 
random (□) phases are also shown for n = 4. The typical spin conhgurations (for 
60x60 lattice) are show in the insets for two different phases. The structured phase 
(ho = 2.0, w = 1.0) and the random phase (ho = 1.0, w = 3.0) are shown for n = 4. 
The conhguration of pinned phase is similar to that for other two cases (uniform 
and bimodal) and could not be shown here due to the limitation of the region of 
space in pinned phase. 
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